FINAL REPORT 

for 

the ADP Program 

The Winds of B Supergiants 

(RITSS contract number 3400) 


by Derek Massa 


This grant has resulted in several publicaitons which were forwarded as attachments to previous reports. In 
addition to these, we have also completed a final publication which has been accepted by Astronomy and Astrophviscs. 
A copy of this final publication is attached. 


1 




Astronomy <V Astrophysics manuscript, no. 
(will bo inserted by hand later) 


Wind variability of B supergiants 

IV. A survey of IUE time-series data of 11 BO to B3 stars 

R.K. Prinja 1 , D. Massa 2 and A.W. Fullerton 3,4 


1 Department of Physics &: Astronomy. University College London, Gower Street, London WCIE 6BT, L .K. 
e-mail: rkpCstar.ucl.ac.uk 

2 Raytheon ITSS, NASA/GSFC, Code 081. Greenbelt, MD 20771, USA 
e-mail: massaCtaotaomona . gsfc.nasa.gov 

3 Department of Physics and Astronomy, University of Victoria, P.O. Box 3055, Victoria, BC V8W 3P6, Canada 
e-mail: awfCpha.jhu.edu 

4 Center for Astrophysical Sciences, Dept, of Physics and Astronomy, Johns Hopkins Universit.y,3400 N. Charles 
St, Baltimore, MD 21218-2686, USA 


Received ; accepted 

Abstract. We present the most suitable data sets available in the International Ultraviolet Explorer ( IUE) archive 
for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 BO to B3 
stars is analysed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously re- 
duced high-resolution spectrograms. The targets include ‘normal’ stars with moderate rotation rates and examples 
of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness 
and range of wind variability and highlights different structures in the winds of these stars. This work emphasises 
the suitability of B supergiants for wind studies, under-pinned by the fact that they exiiibit unsaturated wind 
lines for a wide range of ionization. 

The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in 
individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, 
bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of 
different fundamental parameters, but also different structures may occur in the same star at a given epoch. We 
discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges 
that these phenomena present to theoretical studies of time-dependent outflows in massive stars. 

In addition, SEI line- synthesis modelling of the UV wind lines is used to provide further information about the 
state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product 
of mass-loss rate and ion fraction ( M q , ) is a factor of ~ i.o, when integrated between 0.2 and 0.9 Voo\ it can 
however be several times larger over localised velocity regions. At a given effective temperature the mean relative 
ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low 
velocity regime is discussed in terms of structured outflows. Mean ion fractions are estimated over the B0 to Bl 
spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained 
for the ion fractions of UV resonance lines may reflect the role of clumping in the wind. 
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1. Introduction 

The radiation-pressure-driven winds of O and early-B type 
stars are variable on characteristic time-scales of hours and 
days, as a result of perturbations in the outflow due t.o 
small- and large-scale structure, respectively. Observable 
diagnostics of the former include extended black absorp- 
tion troughs in saturated IJV P Cvgni profiles (Lucy 


1983), X-ray emission (Feldmeier et al. 1997, Kahn et al. 
2001. Cassinelli et al. 2001). and non-thermal radio emis- 
sion (Bieging et al. 1989, Scuderi et al. 1998). The small- 
scale wind structures responsible for these signatures are 
believed to result from a strong instability in the line- 
driven winds, which operates on linear scales smaller than 
the Sobolev length. This instability leads to multiply non- 
monotonic velocity fields in the wind and localised density 
enhancements which are bounded by pairs of forward and 
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reverse shocks (e.g. Owocki et al. 1988; Feldmeier 1 99a). 
barge-scale structure in hot star winds is primarily stud- 
ied via systematic variability of I V resonance and opti- 
cal emission lines. These changes diagnose organised, spa- 
tially extended perturbations in the wind, which may be 
induct'd by magnetic or pulsational processes arising at the 
stellar surface. We document in this paper the nature of 
large-scale wind structure in 1 1 early-type B supergiant 
stars, based on high-resolution International Ultraviolet 
Explorer (IUE) time-series spectroscopy. 

One of the major discoveries made with IUE was the 
systematic nature of variability in the P Cygni-type line 
profiles of OB stars. The key developmental sequence was 
defined in a number of time-series UV investigations (e.g. 
Prinja et al. 1987. Prinja & Howarth 1988. Henrichs et al. 
1988. Kaper 1993, Kaper et al. 1996). A relatively broad, 
shallow absorption feature appears at < 0.3 of the terminal 
velocity (Cco), which then evolves blue- wards over several 
hours and perhaps even days, while narrowing in velocity 
width. The velocity at which these 'discrete absorption 
components' (DACs) are first detected is roughly constant 
for a given star, but may vary from star-to-star. 

Subsequent IU E studies established a connection be- 
tween stellar projected rotation velocity (e e sin(i)) and 
the rate of recurrence and acceleration of DACs (Prinja 
1988: Kaper et al. 1996). More recent IUE monitoring 
campaigns extending over a few weeks have shown that 
the stellar wind lines vary continuously, implying that OB 
star winds are structured over a range of spatial and tem- 
poral scales. They have also shown that variations in the 
profiles modulate on rather long time-scales (i.e. several 
days). Furthermore, the variability evident in £ Per (07.5 
111(f); Henrichs et al. 1998) and C Pup (04 If: Howarth et 
al. 1995) appears to modulate on the stellar rotation rate. 
Rotation is plausibly also the key factor in the repetitive 
wind activity seen in HD 64760 (BO. 5 lb: Fullerton et al. 
1997), HD 93843 (05 III; Prinja et al. 1998) and HD 91969 
(BO la; Massa et al. 1998). Rotationallv modulated wind 
fluctuations are also apparent in late B and early A super- 
giants (e.g. Kaufer et al. 1996), which have much larger 
radii and longer rotation periods. 

The fundamental origin of wind variability on large 
spatial scales in hot stars remains uncertain. The causal 
link between inhomogeneities at. the stellar surface and 
large-scale wind structure has not been directly estab- 
lished observationallv or theoretically. A currently pop- 
ular idea is that the large-scale variations discussed above 
(and some of those discussed in this paper) are spawned 
by photospheric irregularities which cause the wind from 
different longitudinal sectors on the surface to emerge with 
different densities and/or velocities. The surface structure 
may. for example, be due to the action of non-radial pulsa- 
tions or the effects of ordered magnetic fields. Once differ- 
ent adjacent wind streams are initiated, they then meet to 
form co-rotation interaction regions (CIRs). The rotating 
CIRs inevitably evolve into spiral-shaped shock structures 
(Cranmer fe Owocki 1996). In this model, variations in 


(lie wind absorption result from material {lowing through 
a semi-permanent, ro- rot a ting wind structure. 

1.1. Winds of early B-type supergiants 

Ihe stellar winds of early B supergiants (Bis) are par- 
ticularly suitable for a study of wind dynamics and the 
effects of large-scale wind structure since, (i) their UV 
wind lines are often well-developed but unsaturated, (ii) 
t hese diagnostic lines cover a wide range of ionization, 
dose to the dominant stage of the wind (see Table 2 ), 
(iii) it is generally easier to follow the progression of 
large-scale structure in the UV lines of Bis compared 
to O-type stars, (iv) the winds have lower terminal ve- 
locities, so individual structures evolve more slowly and 
the important Siiv AA1400 resonance line doublets are 
mostly decoupled, (v) their strong UV photospheric lines 
are excellent effective temperature-luminosity diagnostics, 
and (vi) the strongest UV photospheric lines sample the 
wind/photosphere 'interface', often displaying wind be- 
haviour. 

Motivated by these attractions, we describe in this 
paper the time-dependent wind characteristics of eleven 
BO to B3 supergiants, based on 16 separate IUE high- 
resolution time-series data sets, involving the inspection of 
638 spectrograms in total. This work extends our knowl- 
edge of the characteristics of variable UV resonance lines 
in line-driven winds to cooler and less luminous stars. The 
data are used to demonstrate the presence of a range of 
empirical wind properties and implied varied perturba- 
tions in the outflows. We report here on the incidence 
and repetition of DACs, evidence for rotational modula- 
tion, ionization-state changes, shock formation, and spon- 
taneous wind enhancements. The diverse wind variability 
signatures in early Bis are compared and contrasted to 
those previously noted in luminous O-type stars. The later 
sections of this paper present SEI line-synthesis modelling, 
which is used to explore further the ionization mixture and 
structured nature of the winds. 

2. The IUE time-series data sets 

Our objective was to carry out a systematic survey of the 
temporal nature of wind activity in early B supergiants. 
We therefore restricted our selection to time-series of high- 
resolution (A/AA ~ 10 4 ) IUE SWP (AA1150 to 1900A) 
data sets. A list of the program stars is given in Table 1 . 
They include luminosity classes la and lb for a spectral 
range BO to B3, and have a wide range of projected ro- 
tation velocities. The corresponding 16 data sets analysed 
are also outlined in Table 1 , and include 4 cases where 
short-term wind variability may be compared over two or 
more epochs separated by a few years. The majority of 
stars in our sample are ‘normal’, slow to moderate rota- 
tors, that are representative of their spectral sub-class. 
Two main exceptions are HD 64760 and HD 157246 (7 
Ara) which were selected to examine the effects of rapid 
rotation on wind structure. 
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'Table 1 . Program stars 


HD 

Sp. Type 

Log ( L + j L,: 

D tf (kK) 

tv sin( i ) 
( km s" 1 ) 

P rot ( max ) 
(days) 

Date of I HE 
time series 

Duration 
( days ) 

Sam pling 
(hours) 

37128 

B0 la X" 

5.78 

28.6 

91 

17.7 

1 980 

0.7 

1 

91969 

BO la 

5.78 

28.0 

83 

14.4 

1996 

30 

9 

L 64402 

BO lb X" 

5.38 

28.6 

* i 

13.3 

1986/93/95 

6/16/16 

3-24/24/21 

204172 

B0 lb X" 

5.38 

28.6 

87 

11.6 

1990/93 

2/16 

3/24 

64760 

B0. 5 lb 

5.10 

23.1 

216 

5.2 

1993/95 

6/15 

3/3 

L67756 

B0. 5 lb 

5.50 

23.1 

79 

14.1 

1986 

5.5 

12 

150168 

BI lab-lb 

5.38 

20.3 

121 

13.3 

1993 

6 

3 

47240 

Bl lb 

4.90 

20.3 

103 

11.3 

1993 

16 

24 

157246 

Bl lb 

4.90 

20.3 

245 

4.8 

1993/95 

6/16 

3/24 

96248 

BC1 la 

5.12 

19.4 

83 

24.3 

1996 

30 

9 

53138 

B3 la 

5.24 

16.3 

58 

45.8 

1987 

13 

12-36 


'/iff and Log (Lm/Lcj) from Humphreys k. McElrov ( 19X4); v e sin (i) from Howarth et al. (1997). Spectral types from Walborn 
(1972, 1976), Garrison et al. (1977), Hiltner et ai. (1969), and Morgan et al. (1955). N - refers to stars whose spectra are 

nitrogen deficient. 


The data were either extracted from the archives main- 
tained at the Goddard Space Flight Center (USA) and 
the Rutherford Appleton Laboratory (UK), or they were 
originally part of our observing campaigns. All the spec- 
tra were extracted in a uniform manner from the two- 
dimensional, photometrically linearised images (PI) using 
the UK Starlink package IUEDR (Giddings et al. 1995). 
One of the advantages of the IUEDR extractions is the 
spectra are highly reliable at the shortest wavelengths, 
thus providing high data quality for the important lines at 
Si ill A1206 and N v AA1240. The wavelength scale calibra- 
tion was improved by centering the spectra on selected in- 
terstellar lines (in an echelle order-dependent correction). 
The individual spectra were binned onto a wavelength grid 
with a regular sampling of 0.1 A. Typically the signal-to- 
noise in the continuum is — 20. Individual wind-formed 
line profiles were normalised to a local continuum by fit- 
ting a straight line through several adjacent ‘line-free’ win- 
dows. 

The mean IU E SWP spectrum for each program star 
is shown in Fig. 1. The principal diagnostic wind and 
photospheric lines available in these spectra are listed in 
Table 2 (Massa 1989). The data sets for the following stars 
have been previously reported in the literature: HD 64760 
(Massa et al. 1995b; Prinja et al. 1995: Fullerton et al. 
1997 ; Howarth et al. 1998), HD 157246 (Prinja et al. 1997), 
and HD 91969 (Massa et al. 1998). These data are pre- 
sented in this paper to place them into wider comparative 
context. 

3. Time-dependent stellar wind phenomena 

The characteristics of variable UV lines in our pro- 
gram stars are presented in this section. Fluctuations are 
present in the stellar wind lines of every B supergiant in 
our sample, and in every time-series data set examined 
here. The variations are primarily described by a gallery 
of grev-scale image representations (dynamic spectra) of 
the time-series data (Figs 2 to 17). This image format 
is an effective method for identifying systematic patterns 


in variable line profiles, and for comparing the evolution 
of these features between different spectral lines. In cases 
where the Doppler shifted absorption profiles are unsatu- 
rated (i.e. the majority in our sample), quotient spectra 
are displayed after being normalised by the mean pro- 
file for the time-series. For saturated or black’ absorption 
troughs in the Civ lines of HDs 37128, 91969, 157246 and 
96248, we display residuals, subtracted from the mean pro- 
file. In all cases darker shades in the images indicate re- 
gions of enhanced absorption (reduced flux) with respect 
to the mean profile. 

3.1. List of observed stellar wind phenomena 

Several empirical phenomena can be identified from the 
images in Figs 2 to 17. They are defined below, prior to 
discussion of their occurrence in individual stars: 

Discrete absorption components ( DACs ) - these are the 
features most commonly associated with variable UV res- 
onance line profiles in luminous OB stars. They usually 
appear as localised (in velocity) absorption enhancements 
that migrate from (at least) intermediate velocities to the 
blue-ward profile edge, with a corresponding decrease in 
velocity width. Typically the accelerations of the DACs 
are < 50% of values for the underlying wind, as deter- 
mined from canonical '(3 ~ V velocity laws. 

Bowed structures - these are optical depth changes that 
appear suddenly over a wide velocity range, and which 
develop simultaneously towards lower and higher blue- 
shifted velocities (see e.g. Fullerton et al. 1997). 

Other features - we use this term to identify either less 
organised ‘bursts’ of absorption that usually last a short 
time only (hours), or more substantial absorptions that 
appear suddenly over a large velocity range. These fea- 
tures do not show any clear evolution in their velocity or 
width as a function of time. 

Near- zero velocity changes - cases where line profile vari- 
ability can be traced to less than 0.05 of the terminal 
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Fig. 1. Mean normalised spectra for the 1U E time-series of the 11 program stars. The important spectral lines are identified in 
Table 2. (An intensity shift of 1.5 has been used in this display between successive spectra.) 


velocity. These fluctuations are due to deep-seated distur- 
bances close to the stellar surface. 

Recurrence - This is simply a descriptive term relating 
to a specific (extended) data set, to draw attention to 


repetitive structures. The repetition is normally due to 
the occurrence of DACs or bowed structures. 

Ionization variation - implies that the ionization state of 
a specific wind feature changes with time. It may for ex- 
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Table 2. Key IV E line diagnostics for early- B supergiants winds 


Ion 

Ao (A) 

l.P. (ev) 

E.P. (eV) 

region 

C ii 

15 5 1.5:5. 15.55.71 

24.4 

0.00. 0.01 

wind 

C IV 

1510.20. 1550.77 

64.5 

0.0. 0.0 

wind 

N v 

1238.82. 1242.80 

97.9 

0.0. 0.0 

wind 

Al m 

1854.72. 18(.»2.79 

28.5 

0.0, 0.0 

wind 

Si ill 

1206.50 

33.5 

0.0 

wind 

Si m 

1296.73 

33.5 

6.55 

wind / photosphere 

Si in 

1312.21 

33.5 

10.28 

photosphere 

Si III 

1417.24 

33.5 

10.28 

photosphere 

Si iv 

1393.70. 1402.77 

45.1 

o 

6 

o' 

d 

wind 


ample become more enhanced in absorption in some lines, 
while exhibiting an increased flux in other spectral lines. 
Ionization stratification - cases where one of the phenom- 
ena identified above are seen at different velocities in dif- 
ferent ions. This is usually an optical depth effect between 
low and high ionization lines. 

Results on the occurrence of these variable wind phe- 
nomena is our program stars are summarised in Table 3. 
The winds of early Bis exhibit diverse patterns of stel- 
lar wind variability. Similar phenomena can occur in stars 
that vary in fundamental parameters, and different wind 
structures can co-exist in the wind of a given star at a par- 
ticular epoch. Some notes on the temporal wind behaviour 
in the individual stars are presented in the remainder of 
this section. 

3.2. Notes on individual stars 

HD 37128 {e Ori, BO la N" ; Fig. 2). A well-studied, mod- 
erately nitrogen weak star (Walborn 1976), part of the 
Ori OBI association, e Ori has been cited as a single-lined 
spectroscopic binary (e.g. Frost 1909; Morrell &: Levato 
1991), however its binarity remains uncertain and, for ex- 
ample, Jarad et al. (1989) found it to be a radial velocity 
constant. Variability in the hydrogen lines has been re- 
ported by Gry et al. (1984) and Ebbets (1982). Previous 
suggestions that it is a variable thermal radio source are 
discounted by Blomme et al. (2002). 

UV variability in the short intensive fUE time-series is 
shown in Fig. 2. We note the progression of a DAC feature 
between — -900 to -1200 km s~ l over — 14 hours. High 
velocity narrow absorption components are also present 
towards the short-ward edge. Strong photospheric lines 
(e.g. C III 1247 and Aim 1855) reveal radial velocity shifts 
of ~ 20 km s” 1 . 

HD 91969 (B0 la; Fig. 3). A member of the Car OBI 
association (e.g. Humphreys 1978). It is one of the bright- 
est members of the NGC 3293 open cluster in the Carina 
region (e.g. Feinstein & Marracu 1980). The IUE time- 
series data studied here were secured during our 4 MEGA 
IF campaign in 1996, and some of our findings have been 
summarised by Massa et al. (1998). This star is an excel- 
lent example of the variety of diagnostic UV lines available 
in early Bis. Massa et al. ( 1998) reported a 7.8-days period 
for the repetition of DAC features propagating blue-ward 


toward the profile edge. These features are seen at very 
low velocities in Aim, and their progression 'connects’ 
with corresponding events seen at higher velocities in Si III, 
Si IV and Nv. There is however evidence for ion stratifi- 
cation between activity in Aim and Siiv. The temporal 
coherency is strongest in the low ion species which form 
in the deeper regions of the wind. The periodic line profile 
changes are seen down to very low velocities (less than 
0.05 i'oq)t together with variations in photospheric lines 
due to Sim, Fe III, and Feiv. The latter may be tracers 
of localised surface inhomogeneities responsible for wind 
activity. 

HD 164402 (B0 lb N“; Figs. 4 to 6). A member of the Sgr 
OBI association (Humphreys 1978), this star has one of 
the lowest projected rotation velocities in our sample. It 
is moderately nitrogen deficient, with no indication of bi- 
nary companionship (e.g. Levato et al. 1988). HD 164402 
provides an opportunity to study systematic wind activ- 
ity over three separate epochs. The 1986 data set (Fig. 
4) has a poorer temporal sampling, but reveals a distinc- 
tive DAC progressing from 1200 km s -1 to ~ —1600 

km s* 1 over 6 days. This feature is accompanied by more 
sporadic events, which are well spread in velocity, but un- 
resolved temporally (e.g. at T ~ 3.5, 4,2 and 5 days, where 
T is time relative to the first observation in the series). The 
1993 time-series (Fig. 5) spans almost 16 days and reveals 
a remarkably slowly migrating DAC, with a mean acceler- 
ation of - 3 x 10 -4 km s~ 2 . Other features with different 
characteristics and more rapid evolution are also noted at 
T ~ 6 and 12 days. Drastic variations are present in the 
1995 data set (Fig. 6) around T ~ 14 days. The high- 
est velocity absorption (at ~ — 1750 km s -1 ) suddenly 
decreases and is accompanied by lower velocity DAC fea- 
tures at ~ —1300 km s -1 and —1000 km s -1 . The sudden 
reversal in optical depth across the lines may represent 
the signature of a large scale shock. 

HD 204172 (69 Cyg; B0 lb N“; Figs. 7,8). A moder- 
ately nitrogen weak star (Walborn 1976), belonging to 
the Cyg OB2 association. Lennon et al. (1992) highlight 
several discrepancies between the B0 lb spectral type as- 
signment and the appearance of the optical spectral line 
profiles, including H lines being too narrow and Siiv too 
strong. They suggest a revision to BO. 2 la. The UV reso- 
nance lines of 69 Cyg are generally stronger than those of 
HD 164402, thus favouring a slightly more luminous spec- 
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Table 3. Observed stellar wind phenomena 


HD 

DAC 

Bowed 

structure 

Other 
feat ure 

near-zero 

var. 

ion 

strat. 

ion 

var. 

recurrence 

3712s 

V 







9 1 969 

V 


/ 

V 

V 

i 

yj 

V 

10-1-102 

v" 


V 

■* 



- 

204172 

■* 


%/ 





64760 

V 


sj 

V 

V 

yj 

v' 

167756 

V 


yj 

V 

v 7 


- 

150168 

V 


c 

V 

yj 

yj 

y/ 

47240 

V 


yj 

V 

V 

- 

- 

157246 

V 


c 

- 

yj 

yj 


96248 

V 

- 

y/ 

- 

yj 

- 

- 

53138 

sj 


yj 

- 

- 

- 

- 


tral type. The densely sampled 1990 time series (Fig. 7) 
reveals a gradual increase in profile absorption strength 
over ~ 2.2. days. The fluctuations extend over a broad ve- 
locity range that cover almost the entire Si IV absorption 
trough. Corresponding changes are seen at much lower ve- 
locity (extending close to zero velocity) in C II A 1335 and 
A 1 III A 1855. The velocity differences between low and high 
ions are principally an optical depth effect. There are also 
traces of DAC features at ^ —1700 km s -i (e.g. T ~ 0 to 
0.5 days) in Siiv, Civ and N v. Some changes are appar- 
ent in the strong Cm A1247 'photospheric' line. The 1993 
data set (Fig. 8) is more extended and covers almost 16 
days. Two enhanced absorptions are present over a sub- 
stantial velocity range. These may be similar events to the 
phenomena seen with greater time resolution in 1990 (Fig. 
7). There is a tentative indication that the feature recurs 
over ~ 12 days, which is close to the estimated maximum 
rotation period of the star (Table 1). Some variations are 
also noted in the deep-seated Cm A1247 and Sim A1300 
triplets. 

HD 6^760 (B0. 5 lb; Figs. 9. 10). This is one of two rapid 
rotators in our sample. The data sets for the two epochs 
shown here were discussed in detail by Massa et al. (1995b) 
and Fullerton et al. (1997). HD 64760 exhibits a variety of 
spectral features that indicate co-existing wind structure. 
In fact, every category of stellar wind phenomena listed in 
Table 3 is present in this star! The ‘MEGA I' campaign 
during 1995 (Fig. 10) in particular reveals a striking 1.2- 
day modulation, where individual features cover a large 
velocity range, and have a bow-shaped morphology such 
that they extend blue-wards and red-wards at the same 
time. Ionization variability is indicated by the fact that the 
N v data appear to have nearly twice the modulation fre- 
quency of Si IV (Fullerton et al., 1997). Slower migrating, 
high velocity DACs are also present, which persist over 
several days. The recurrence time-scale for the DACs is 
not constrained. An exceptionally strong absorption event 
is present in high-ion species in the (shorter) time series 
data for 1993 (Fig. 9, around T — 2 days). It is accom- 
panied by variations essentially down to zero velocity in 
e.g. Cll AA1335 and Al ill A 1855 * plus disturbances in the 
*photospheric’ Si m A 1300 triplets. The onset of these sub- 


stantial changes is also marked by changes in the relative 
ionization mixture, such that the overall ionization of the 
wind increases during this event (Massa et al. 1995b). It 
is possible that this strong event, and the overall enhance- 
ment seen at T — 4 days (Fig. 9), are similar to the bowed 
structures apparent in the more extensive MEGA I data 
set (Fig. 10). 

HD 167756 (B0. 5 lb; Fig. 11). A distant halo star, located 
about 4 Kpc away at a Galactic altitude of —0.85 Kpc 
(e.g. Savage & Massa 1985). Ionization shifts are appar- 
ent in the wind-formed lines, with the high velocity (^ 

— 1500 km s" 1 ) DAC seen in Siiv (T ~ 0 to 3 days) only 
being apparent at less than —500 km s -1 in Aim. The ab- 
sorption enhancement seen in Siiv, N v, and Civ towards 
the end of the observing run (T > 4 days) may also be 
accompanied by variations close to zero velocity in Si III 
A1207 and Aim A1855. 

HD 150168 (B1 Iab-Ib; Fig. 12). A member of the Ara 
OBI association (Humphreys 1978). The UV time-series 
reveals a radial velocity motion in strong UV photospheric 
lines (e.g. Cm A1247) due to binary motion. Several in- 
teresting wind effects are apparent in this star. Near the 
start of the time-series there is a maximum in absorption 
in N v which corresponds to a minimum in Siiv (i.e. ion- 
ization variability). Close to T ~ 4 days a bowed structure 
is present, with the absorption minimum in N v lagging 
the Siiv minimum (i.e. ion stratification). A much weaker 
bowed feature may be present at T ~ 2 days. Multiple 
high velocity DACs (> —1500 km s _1 ) are seen in Siiv, 
C iv and N v. 

HD 47240 (B1 lb; Fig. 13). A member of the Mon OB2 
association. Binary radial velocity motion is seen in strong 
UV photospheric lines (e.g. Cm A1247).The double DAC 
feature seen at the start of the time-series at ~ —700 km 
s"* 1 in Siiv is present at a lower velocity (< —500 km s^ 1 ) 
in C II and Al III. This enhancement is seen at almost zero 
velocity in Civ, however there is very little line flux in 
Si ill A1207. The variations in Siiv between ~ —500 and 

— 1000 km s” 1 may be repetitive over a 3 to 5 day time 
scale. 

HD 157246 (7 Ara, B1 lb; Figs. 14,15). This star has 
the highest projected rotation velocity in our sample of 
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stars (Table L). The 1993 UV time-series was discussed by 
Prinja et al. (1997: see Fig. 14). The data reveal a com- 
plex pattern of co-existing discrete features at low veloci- 
ties (red-ward of —750 km s -1 ) and high velocities (blue- 
ward of -1500 km s" 1 ). These regions are ‘divided' by 
the appearance of a very sharp (in velocity width) DAC- 
like feature. The location of this narrow feature is shifted 
by ~ 400 km s" 1 between Siiv and Aim (or (Til). The 
two-component wind properties, and overall profile mor- 
phologies, suggest evidence for an equatorially compressed 
wind, with structure in the polar components (Prinja et al. 
1997). The narrow DAC-like feature is not seen in the 1995 
IUE time-series (Fig. 15), and the high speed wind be- 
tween -1000 and -1500 km s _1 is also no longer present. 
The line profiles in the 1995 data are more representa- 
tive of the normal state of 7 Ara (e.g. when compared to 
isolated archival spectra taken at earlier epochs). Some 
ionization state variability is observed between Si IV and 
N v (see e.g. the fluctuations between T ~ 2 to 8 days in 
Fig. 15). In both data sets variations are present down to 
near-zero velocities in the low ions. 

HD 96248 (BC1 la; Fig. 16). Carbon enhancement has 
been reported in this star by Dufton (1972) and Walborn 
(1976). HD 96248 is a radial velocity constant (e.g. Levato 
et al. 1988). The ~ 30 day time-series shown in Fig. 16 is 
from the 'MEGA IT Campaign of 1996. The Si IV and C IV 
profiles are nearly saturated, and the line profile variability 
is confined to the blue edges (i.e. > —750 to —1000 km 
s _1 ). A gradual increase in absorption in the blue edge 
occurs from the start of the time-series to almost the end 
of the run, only weakening slightly after T — 24 days. 
These variations are mimicked in the low ions (Aim and 
Cll), but at slightly lower velocities. Additional variations 
due to DACs are also seen between — —200 and —600 km 
s" 1 in Aim, but they do not evolve into the structure seen 
at high velocities. The blue edge in Cll A1335 is not as 
shallow or extended as in the other lines. 

HD 53138 (o 2 CMa, B3 la; Fig. 17). A member of the Coll 
121 association. HD 53138 exhibited very complex and 
highly variable characteristics in the Ha profile survey 
of Ebbets (1982), with variable absorption and emission 
components. Low velocity (red-ward of —500 km s -1 ) UV 
variability in present in the ~ 13 day time series shown 
in Fig. 17. A DAC feature is present between T ~ 0 to 
4 days, which is essentially confined to high velocities. 
Various sporadic changes are apparent later in the run. 
Some degree of ion stratification is present, such that Cll 
variations during the first 4 days appear at lower velocities 
than in Al III. 

4. Comparisons 

It is clear from Figs. 2 to 17 (and Table 2) that the stellar 
wind activity in early B supergiants, as revealed by the 
UV resonance lines, is substantial and has highly varied 
characteristics. 

It is interesting to compare the BI wind activity pre- 
sented here with that of 10 O-type stars described by 


Kaper et al. (1996), also based on high-resolution IUE 
spectroscopy. The O stars generally do not exhibit such 
a variety of temporal phenomena, and in particular show 
limited evidence for ionization related effects. (The com- 
parisons are restricted, however, since Si IV is often the 
only reliable indicator of wind activity in 0-tvpe stars, 
as Nv and Civ are saturated and the sub-ordinate N iv 
1718 is only weakly present.) The fundamental character- 
istic patterns of variability in the O star data are in the 
form of blue-ward migrating DACs, though their proper- 
ties may differ a lot between e.g. acceleration rates, mul- 
tiplicity, and line-of-sight velocity dispersions. The coher- 
ence of systematic structure is generally greater in the 
UV lines of Bis compared to those of 0 stars. In terms of 
clearly defined - obviously modulated - wind structure, 
the B supergiant examples of HD 64760 and HD 91969 
described above, are best matched by the 07.5 III stars £ 
Per (HD 24912) and 68 Cygni (HD 203064). 

In order to quantify the extent of variability as a func- 
tion of velocity, we computed the temporal variance spec- 
trum (TVS, see Fullerton et al. 1996) for the Si iv line 
profile for all our program stars: 


(TVS),- = 




(i) 


where 5, 7 is the normalized intensity of the ith pixel 
in the jth spectrum, S,- is the weighted mean of the 
normalized intensity, aj c is the inverse of S/N of spec- 
trum j measured in an adjacent continuum band, and 
t — 1 

. The results are shown as root 


I* 


CJ 


mean square percentages (= TVS^ x 100) in Fig. 18. In 
all the stars with extended time-series data (i.e. spanning 
more that a few days), the absorption variability in Si iv is 
observed to extend down to less than 0.1 t>co, and to essen- 
tially rest velocity in a few cases (e.g. HDs 64760, 150168, 
47240 and 157246). This is a lower velocity than is gener- 
ally the case for O-type stars, and the result suggests that 
the winds of B stars are structured differently. Greater 
variability in the lower velocity portion of the P Cygni 
profiles in Bis indicates a lack of forward-scattered emis- 
sion compared to O stars. This could be due to a mul- 
tiply non-monotonic velocity law (which implies a higher 
degree of structure) or large-scale deviations from spheri- 
cal symmetry (which decrease the net amount of forward 
scattering) or clumping (less scattering in all directions). 

The maximum amplitude of the variations in Fig. 18 
mostly occurs between ~ 0.7 to 0.9 Uoo- HD 91969 and 
HD96248 exhibit the largest amplitude changes in Siiv. 
Though the variance diagnostics in Fig. 18 reveal no case 
of substantial variation in the emission components of the 
P Cygni profiles, weak systematic changes have been noted 
in the UV emission lobes of HD 64760 (e.g. Fullerton et 
al. 1997). 

Measurements were also made of the ratios of mini- 
mum flux to maximum flux as a function of velocity for 
selected lines in the time-series data of each star. The over- 
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Fig. 2. Grey-scale representation of variability in HD 37128. 



Fig. 3. Grev-scale representation of variability in HD 91969. 
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Fig. 4. Grey-scale representation of variability in HD 164402 (1986). 
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Fig. 5. Grey-scale representation of variability in HD 164402 (1993). 
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Fig. 7. Grey-scale representation of variability in HD 204172 (1990). 



R.K. Prinja, D. Massa anti AAV. Fullerton: Wind variability of B supergiants 


11 


HD 204172 - BO lb - 1993 


Si IV A 1400 



0 L_jl 


-2000 0 2000 


C IV A 1550 N V A 1240 



-2000 0 2000 -2000 0 2000 


Si m A 1207 




A1 m A 1865 



I L_ 

-2000 0 2000 


Velocity (km/s) 


Fig. 8. Grey-scale representation of variability in HD 204172 (1993). 
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Fig. 9. Grey-scale representation of variability in HD 64760 (1993). 
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v-scale representation of variability in HD 64760 (1995). 
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Fig. 11. Grey-scale representation of variability in HD 167756. 
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Fig. 12. Grey-scale representation of variability in HD 150168. 
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Fig. 13. Grey-scale representation of variability in HD 47240. 
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Fig. 14. Grey-scale representation of variability in HD 157246 (1993). 
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Fig. 15. Grev-scale representation of variability in HD 157246 (1995). 
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Fig. 16. Grey-scale representation of variability in HD 96248. 
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Fig. 17. Grev-scale representation of variability in HD 53138. 
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Fig. 18. Temporal variance (RMS) for the Si IV spectral lines in our program stars. The velocity scale is with respect to the red 
component of the doublet, and the rest wavelength of the blue component is marked by an arrow. The horizontal dashed lines 
mark the level of significance for a probability of 95%. 


all trend in N V\ Siiv and Civ is of a gradual decrease in 
Fmin(t } )/F m ax(t’) (i.e. greater relative absorption or opti- 
cal depth) with increasing velocity. The mean ratios (over 
all stars) vary as: F min (t J )/Fmax(^)[N v] — 0.81 (at 0.2 tv^) 


to 0.68 (at 0.7 Voo ), F min (v)/F max (t/)[Siiv] - 0.80 (0.2 
t? co) to 0.33 (0.7 Voo), and F min (i;)/F max (v)[Clv] - 0.53 
(0.2 t’oo) to 0.05 (0.7 Voq). The most extreme variations in 
lines other than Siiv (Fig. 18) occur in HD 150168 and 
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HD 47240 in N v and Civ. The weaker A! Ill flux ratios do 
not exhibit a clear trend in velocity and mean values vary 
between — 0.75 to 0.9 - stronger absorption variations are 
seen in this line in HDs 47240. 15724(5 and 9(5248. 

4.1. Rotational modulation 

Aside from the data sets of HD 91969 (Fig. 5) and 
HD 64760 (1995; Fig. 10), the remaining IU E time-series 
data of the program stars do not reveal unambiguous evi- 
dence for repeatable or cyclic modulations. In several cases 
this is likely because the IU E observations are not ex- 
tended enough to cover stellar rotational time-scales. We 
nevertheless performed a Fourier period-finding analysis, 
using the CLEAN algorithm (Roberts et al. 1997) to itera- 
tively remove the window function from the Fourier spec- 
trum. We examined all the N V, Si iv and Civ time-series 
data sets for the program stars (where the line profiles had 
been previously normalised to the local continuum). The 
results indicate that (excluding HD 91969 and HD 64760) 
the only other cases of marginal significance where power 
is apparent across a reasonable fraction of the absorption 
trough (i.e. over few hundred km s~ J ) are (i) a modula- 
tion of ~ 4 days in the 1993 data set of HD 204172 (Fig. 
8), and (ii) fluctuations on ~ 1 day in the N v data set of 
HD 150168 (Fig. 12). 

Variability in the Si IV lines is compared more directly 
in the images shown in Fig. 19, where the stars are or- 
dered in terms of estimated maximum rotation period. 
There is perhaps some overall similarity between the mod- 
erate rotators HDs 47240, 204172, 164402 and 91969. Only 
HD 150168 shows bowed features that resemble the strik- 
ing patterns seen in HD 64760. There is also a resem- 
blance between the very slowly propagating, high veloc- 
ity DAC seen in HD 64760 (blue-wards of — -1200 km 
s -1 ) and DACs seen towards the short-ward profile edge 
in HD 204172 and HD 164402. 

In O-type stars the UV resonance line profile variabil- 
ity is usually dominated by the biue-ward migrating, re- 
current, DACs. The apparent connection between modu- 
lated wind activity and stellar rotation in these stars is 
then based on measured trends in the recurrence and evo- 
lution time-scales of the DACs. We have show n in Sections 
3 and 4 that the line profile behaviour of Bis is more com- 
plex, such that the variance (e.g. Fig. 18) is generally not 
due to a single phenomenon such as DACs. but instead 
reflects the contributions of different variable structures. 
Furthermore, the empirical properties of these co-existing 
features do not all follow the same evolution time-scales. 
Examples of fairly localised (i.e. velocity widths ~ 0.1 to 
0.3 t’co) blue-ward migrating features that resemble DACs 
in O-type stars are seen in our sample in HDs 91969, 
164402, 64760 (1995), 150168, 157246 (1993) and 96248. 
The measured mean accelerations for the DACs in the 
highest v e sin (i) stars, HDs 64760 and 157246, are — 90 
km s' 1 d' 1 and 140 km s -1 d -1 , respectively. The corre- 
sponding values for the slower rotators, HDs 96248, 150168 


and 1(54402 are — 7 km s _1 d -1 , 35 km s _l cl -1 and 20 
km s* 1 d“ l . Therefore, aside from the other variability 
phenomena identified in Section 3.1. these results are con- 
sistent with the view that the time-scale of stellar wind 
activity in radiation- pressure-driven, luminous earlv-type 
stars is related to the rotation period of the underlying 
star. However, our study also highlights other broad ab- 
sorption features which can arise and evolve substantially 
over comparatively short time-scales (~ 1 to 2 days) in 
seemingly moderate or slow' rotators, such as HD 204172 
and HD 150168 (see e.g. Fig. 19). 

5. Profile fitting with the SEI code 

In this section we present the results from wind line pro- 
file modelling, to derive further details of the UV vari- 
ability characteristics and ionization states of early B su- 
pergiants. Specifically, our objective was to derive values 
for the product of mass-loss rate and ion fraction (M <?, ), 
for the mean line profiles of the time-series and for cases 
of overall maximum and minimum observed absorption. 
We also examined trends in the relative ionization ratios. 
Our approach in this analysis is limited to the simplified 
method of employing steady-state line synthesis models 
to match observed profiles that are obviously affected by 
temporal wind structure and inhomogeneities in the out- 
flow. Nevertheless, some useful information can still be de- 
rived from the success and failures in matching line profiles 
due to different ions. 

We modelled the N v, C iv, Si IV and Si III (A1207) wind 
lines (and in a few cases also Aim and Cll) using the 
'Sobolev with exact integration’ (SEI) method described 
by Lamers et al. (1987). In addition, we adopted the modi- 
fied treatment of the wind optical depth law introduced by 
Massa et al. (1995b), such that the radial optical depth 
Trad(^) is treated in 10 independent velocity bins (each 
0.1 Voq wide). This approach provides greater flexibility in 
matching absorption optical depths in a wind that is af- 
fected by time-variable structure. We assume a standard 
parameterised law for the wind of the form 

w = wo 4- (1 - u’o)(l - 1/*)^ (2) 

where w = v/v <*, and x = r/R, where R is the stellar ra- 
dius. We select 0 = 1.0 and wo = 0.01, and generally find 
that the model profiles do not differ grossly for values of 
d between 0.5 to 2.0. A range of 0.05 to 0.1 was adopted 
for i?turb i which is the small-scale velocity dispersion pa- 
rameter (i.e. the localised velocity dispersion is typically 
~ 100 km s“ 1 ). 

Since the UV photospheric spectrum of early Bis can 
be fairly strong, it is necessary to allow for their effects 
in the SEI profile fitting. The high-resolution IU E pho- 
tospheric spectra of the following stars were adopted (by 
spectral type) as ‘standards* and explicitly included as 
a lower boundary for the calculations: HD 36512 (B0 
V), HD 36960 (B0. 5 V), HD 31726 (B1 V), HD 36959 
(B1.5 V), HD 35468 (B2 III), HD 207330 (B2.5 III) and 
HD 209008 (B3 III). These ‘standard’ stars were selected 
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Fig. 19. Comparison of the Si iv variability characteristics of selected program stars; the estimated maximum rotation period 
is increasing from left to right. 


as low v e sin (i) targets, covering the range of spectral 
types of our program stars, and exhibiting very weak or no 
evidence for winds in their resonance lines. The input pho- 
tospheric spectra were convolved with a rotational profile 
to match the rotational broadening of the program stars, 
and generally the agreement between the corresponding 
photospheric absorption lines was very good. 

With the parameters in the velocity law and photo- 
spheric input fixed, the r ra d(m) bins were the only free pa- 
rameters in the line fitting scheme. The procedure adopted 
was to fit the individual profile, bin by bin, progressing 
from the blue-ward profile edge to less negative veloci- 
ties (see Massa et al. 1995b). We modelled in this manner 
the mean Nv, Civ, Siiv, and Sim profiles for the time- 
series data of each program star. In addition, to gauge 
the range in parameters due to profile variability, we also 
modelled for Nv, Si iv, and C IV, the representative (over- 
all) maximum and minimum absorption spectra (the IUE 
SWP images used in these cases are listed in Table 4). 
The Si III A1207 lines are mostly too weak and their pho- 
tospheric absorption lines are less well matched by the 
standard stars, to determine variability characteristics by 
profile fitting. The mass-loss rate results are summarised 
in Table 4, where the M qi values are integrated over 0.2 


to 0.9 i^oo . Examples of the fits to the mean profiles are 
shown in Fig. 20. We note the following key points from 
this analysis. 

The overall quality of the line profile matches is 
good, which partly reflects the flexible optical depth ‘law 1 
adopted. It is notable that mostly in Civ - and to a 
lesser extent in Si IV - the low velocity absorption (i.e. 
below about 0.25^) is not well matched. The spherically 
symmetric models often predict excess forward-scattered 
emission in the low-velocity part of the P Cygni absorp- 
tion trough. The discrepancy scales qualitatively with the 
strength of a given spectral line, such that it is gener- 
ally worse for the maximum absorption cases than the 
minimum ones for a single star. This does not hold as a 
function of mass-loss rate though, and for example the 
Si iv mismatch is just as large for HD 150168 as it is for 
HD 37128, for which M qi is a factor of ~ 5 larger (Fig. 20). 
Puis et al. (1993) discuss in detail the fact that reduced 
emission results when account is taken of back-scattering 
in several different resonance zones in structured winds. 
The reduced emission acts so as to yield an apparent in- 
creased absorption in the low velocity region between ~ 0 
and 0.3t;<x>. Also if the winds are not spherically symmet- 
ric, then the models will produce too much scattered light 
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Table 4. Results from SK! model profile fits 




A/ 7(<V+ ) 

Mq(Si' + ) 

M 'i(C t+ ) 

\f q{Sr + ) 

Star | HD) 

i’x 

mean 

range 

mean 

range 

mean 

range 

mean 

37128 

1700 

2.75(-9) 

2.60-3. 10( -9) 

1.64(-8) 

1. 13 - 1.79(-8) 

>4. 791-9) 


1 . 56( -9 ) 

91969 

1500 

4.31 (-9) 

4.00-5. 25(0) 

1.60(-8) 

0.88-1.60(-8) 

> 4 . 0 1 ( - 9 ) 


1.39(-9) 

1 G-H 02 

1650 

6.3K-L0) 

5.2l-6.54(- 10) 

l.UI-0) 

l. L2-1.66(-9) 

1.07(-9) 

0.82 -1.30 (-9) 

1.34(-10) 

204172 

1750 

8.36(- 10) 

0.74-1. 53(0) 

3.06(-9) 

2.14-4.34(-9) 

1.92 (-9) 

1.84-3.48(-9) 

3.12(-10) 

64760 

1300 

9.90(- 10) 

0.86-1. 45(0) 

2.88(-9) 

1 . 77—4 . SO ( -9 ) 

1.82(-9) 

1.54-2.41 (-9) 

2.22(-10) 

167756 

1750 

1.52(-9) 

1.46-1. 64(-9) 

3. 25 (-9) 

2. 4 2-3. 37 (-9) 

2.01 (-9) 

1 .93— "2. 2 1 ( -9) 

3. 70(- 10) 

150168 

1370 

5.30(-9) 

4.7G-6.98(-9) 

3.58(-9) 

0.80-5. 76(-9) 

2.34(A) ) 

1.55-2.80(-9) 

2.77(- 10) 

47240 

1000 

2.67(10) 

2.65-3.87(- 10) 

2.57(-9) 

2.20-5. 43(-9) 

> 1.28(-9) 

- 

6.49(-10) 

157246 

600 

7.30( - 10) 

7.00- 14.30(- 10) 

9.801-10) 

5.70-16.50(-'10) 

>4.61(- 10) 

- 

- 

96248 

800 

- 

- 

l.46(-9) 

1.45-1.46(-9) 

5.81 (-9) 

4.29-6.03(-9) 

- 

53138 

520 

- 

- 

2.78(-9) 

2.55-3.07(-9) 

7.27(-10) 

7.20-7.97(-10) 

- 


The products of mass-loss rate and ionization fraction ( M q x ) are derived from fits to the mean profiles. (X.XX(-V) = X.XX x 
10”^ A/iti yr“ l ). The range in this value determined from cases of overall minimum and maximum absorption is also listed. 
The IUE SWP images corresponding to the cases of minimum and maximum absorption that were modelled are: HD 37128 
(8100. 8130), HD 91969 (57187, 57426). HD 164402 (54200, 54225), HD 204172 (38957, 48914), HD 64760 (53522, 53709), 
HD 167756 (29001, 29049), HD 150168 (47098, 47198), HD 47240 (48953, 48930), HD 157246 (54172, 54103), HD 96248 
(57069, 57204), HD 53138 (30153, 30264). We also obtained a few results from fits to Aim - HD 204172. 2.72(-9); HD 47240, 

2 . 4 7 ( -9 ) ; HD 96248, 3.38(-8) and Cii - HD 53138, 6.84(-10). 


at low velocities. Note, however, that the discrepancy at 
low velocities in a fast rotator like HD 64760 is no greater 
than that for an apparent moderate rotator of the same 
spectral type, HD 167756 (Fig. 20). It is also possible that 
to some degree the mis-match at low velocities reflects 
clumping in a highly structured wind so that parts of the 
stellar disk are being obscured by very different amounts 
of material than others. The low velocity mis-match is 
generally much lower for N v compared to Si IV and Civ 
(Fig. 20), implying then that the region occupied by the 
N v emission may be more uniformly filled. 

The mean values of M qi listed in Table 4 are consis- 
tent with the dependence on luminosity for B supergiants 
reported by Prinja fc Massa (1997), who presented prelim- 
inary results for a wider, time-averaged survey. Typically 
the range in M qi implied by the variability between min- 
imum and maximum absorption line profiles is a factor of 
~ 1.5. A clear exception is the factor of 7 change noted 
in M q(Si 3+ ) for HD 150168; this star showed strong evi- 
dence for ionization variability in the time-series data (e.g. 
Fig. 12 and see below). This quantity describes the change 
in structures integrated over a velocity range from 0.2 to 
0.9 too, and is therefore intended to reflect variations in 
the global mass-loss rate. However, the localised variations 
(i.e., in a given velocity bin) can be several times larger. 
For example, at 0.7 Vco (cf. temporal variance in Fig. 18), 
the product of M q(Si 3 +) may vary between a factor of 
4 to 15 for supergiants in the spectral range B0 to BO. 5. 
These changes reflect the degree of structure in the wind. 

Our sample of stars is too small to derive useful infor- 
mation on the temperature dependence of the relative ion 
ratios. However, at a given effective temperature the val- 
ues of q(Si 3 + /q(N 4 +), or q{Si 3 + )/?(C 3+ ), can differ by up 
to a factor of 5. We find that there are no significant differ- 


ences (or trends) in an individual star between the relative 
ion fraction derived from the mean profiles, minimum ab- 
sorption profiles and maximum absorption profiles. Over 
the spectral range BO to B 1 (la and lb) the mean ion 
ratios are q{Si 3+ )/q(N 4 +) ~ 2.5, q(Si 3 + )/q{C 3+ ) ~ 1.7 
and q{Si 2 +)/q(Si 3 +) ~ 0.1. HD 150168 and HD 157246 
(7 Ara) exhibit the most exceptional mean ion ratios of 
q{N 4 +)/q(Si 3+ ). In HD 150168 the ratio is almost 3 times 
larger than for other stars of similar parameters. The dif- 
ference is linked to the ionization change (of unknown 
origin) that occurred during the run (see Sect. 3.2); the 
increase in the ion ratio points to a temporal ionization 
shift from a 'normal state' to a 'high state'. The ratio of 
q(N 4+ )/q{Si 3+ ) in 7 Ara is also several times larger than 
for other Bis. This result may reflect a nitrogen rich stellar 
surface due rotational mixing (in a rapid rotator). As ex- 
pected, the q(C 3 +)/q(Si 3 +) ratio of the BC star HD 96428 
is abnormal, and the mean ratio is between 4 to 8 times 
larger than the other stars in our sample. 

The behaviour of the ionic ratios as a function of ve- 
locity is shown in Fig. 21 for our program stars (based 
on fits to the mean profiles of the time-series). The 
q(Si 2 +)/q(Si 3 +) ratio does not depend on abundance and 
its observed trend indicates a decrease in ionization state 
of the wind as a function of velocity. The relative ion mix- 
ture of q (Si 3+ )/q(N 4+ ) seems to be more stable between 
~ 0.4 to 0.8 Uoo (or ~ 2 to 5 stellar radii for the adopted 
empirical velocity law r ). 

HD 37128 (e Ori) is the only star in our sample with 
a measured radio flux, which may be used to provide 
a reliable estimate of the mass-loss rate. Blomme et al. 
( 2002 ) have conducted a radio and sub-mm study of e 
Ori, and they derive a value of M = 1.9 x 10~ b M yr“ l . 
The (mean) results in Table 4 therefore give ion frac- 
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167756 


150168 
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-10 12 -10 12 

Normalised Velocity (v/v„) 


Fig. 20. SEI model line profile fits for N V, Civ. Si IV, and Si III lines in our program stars. Examples are shown for the cases 
of constructed mean profiles (see also Table 4). 

tions for this star of <p(jV 4+ ) ~ 1.5 x 10~ 3 , g(S* 3+ ) — to Bi of q(N 4 * ) — 2 x 10 “ 3 , </(Si 3 +) ~ 5 x 10~ 3 , q(C 3+ ) 
9.0 x 10 -3 , (/(C 3 " 1- ) > 2.6 x 10~ 3 . and q(Si~+) ^ 8.0 x ~ 3 x 10“ 3 . and q{Si 2 +) ~ 1 x 10 -3 . The line-fitting 
10“ 4 , Alternatively, using the theoretical mass-loss recipe analysis suggests therefore that none of the spectral lines 
of Vink et al. (2001) to determine the mass-loss rates for modelled here represent the dominant stage of ionization 
our target stars, yields mean ionization fractions over BO for early B supergiants (over the UV line-formation re- 
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gion). This is particularly surprising for the case of Si 34 ' 
and C 3+ in the earliest B stars studied, such as e Ori. As in 
0-tvpe stars, the other likely candidates for t he dominant 
ion stages are Si 4+ . C i+ , and N' i+ . 

It may be that the low ionization fractions derived here 
are due to the effects of clumping in the stellar winds 
(e.g. Massa et al. 2002, in preparation). If the winds are 
principally composed of optically thick clumps separated 
by transparent 'voids’, then each constant velocity sur- 
face will only partially contribute to the emission line. 
Furthermore, the face of the star would not be completely 
covered by optically thick material at any velocity. This 
means the absorption part of the profile will not be very 
deep, thus also giving the impression of an unsaturated 
line. As a result, when fitting the spectral line with a ho- 
mogeneous wind model, vve derive a finite optical depth 
(and lower values of M qi ), even though the optical depth 
of the wind might be much larger if the wind material was 
homogeneously distributed. 

6. Discussion 

We have shown that the stellar wind activity in early B 
supergiants is rich and varied. Figures 2 to 17 illustrate 
the fact that a selection of temporal signatures can occur 
in an individual star, and between stars. Different struc- 
tures can also occur in the same star at different epochs, 
though 7 Ara (HD 157246) is the only extreme example of 
this behaviour in our sample. The physical connection be- 
tween these perturbations remains unresolved. There are 
two general forms of variability characteristics that appear 
to co-exist in the UV data in some stars : 

1. Firstly there are broad absorption features which have 
a large spatial coherence (assuming monotonic veloc- 
ity laws) and are relatively short-lived. Examples of 
these can be seen in HDs 64760 and 150168 as bowed 
features, and HDs 91969, 164402. 204172 and 47240 as 
blue-ward accelerating absorption enhancements. 

2. In most of the stars listed above the broad absorption 
features are sometimes ‘accompanied’ by DACs, which 
are characteristically different in that they have longer 
temporal coherence. These features are much nar- 
rower in velocity, accelerate blue-wards more slowly, 
and are quickly confined to the highest velocities near 
the short-ward profile edge. Their origin is enigmatic 
partly because the UV time-series are mostly not ex- 
tended enough to monitor the recurrence time-scales 
(and properties) of these features. It is therefore diffi- 
cult to establish (or rule out) their connection to the 
extended absorption features (or bowed structures). 

The first of the two phenomena are usually linked with 
the effects of CIRs (e.g. Owocki et al. 1995; Cranmer 
and Owocki 1996) passing through the line-of-sight. 
Specifically, localised features or inhomogeneities at the 
stellar surface are predicted to alter the radiative force, 
which leads to the formation of high-density low-speed 
gas streams. These streams are then rammed by the faster 




Fig. 21. The relative ion ratios (qt/q>) are plotted as a func- 
tion of velocity, based on SEI model fits to the mean line pro- 
files of our program stars. 


moving ‘underlying’ wind as the star rotates, leading to a 
spiral pattern. The enhanced absorption identified above 
arises from an extended layer of near-uniform velocity that 
marks the initial response of the unperturbed wind to the 
slower CIR ahead. 

A further possibility then is that the slowly migrating 
DACs (item 2 above) are also related to CIRs, but orig- 
inate in the radiative-acoustic kinks which are predicted 
to trail the CIRs in the hydrodynamical simulations of 
Cranmer & Owocki (1996). 

It is important to realise that both observational phe- 
nomena highlighted above are not always present in every 
OB star. In many O-type stars for instance, only the long- 
lived DAC features are present, i.e. without evidence for 
co-existing ‘bowed’ structures or periodic absorption mod- 
ulations. In this respect, we believe the model developed 
by Hamann et al. (2001) is more appropriate for a star 
with a distinctly structured wind such as HD 64760 (Fig. 
10) or £ Per (Kaper et al. 1997), than for £ Pup where 
the dominant (19.2-hr) variability pattern is due solely to 
blue-ward migrating DACs. 

The strong ionization changes seen in HDs 150168, 
91969, 64760 and 157246 (Table 3) are likely shock sig- 
natures, though we cannot easily assign them to the same 
physical origin. If the interface between slow’ and fast wind 
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streams in the ( ’IR model is associated with strong shocks, 
then we can expect some ionization changes in the I V 
lines as the structure crosses the line-of-sight to the stel- 
lar disk. In cases where the variability phenomenon is seen 
at different velocities in different ions (Table 3), this be- 
haviour may also relate to t he action of spatially confined 
shock regions in the wind, where the post-shock gas is 
denser, lower speed, and less highly ionized than the pre- 
shock material. Multan (1984) and Cranmer Ac Oworki 
(199(5) discuss scenarios which may potentially result in 
velocity jumps of several hundred km s~ 1 . 

An important result from our study is the variety of 
wind activity patterns that are present in the U V time- 
series data of Bis. However, this diversity of behaviour 
does not necessarily indicate that there is a correspond- 
ing variety of fundamental processes responsible for gen- 
erating wind structure. The combination of stellar rota- 
tion rate and viewing (aspect) angle to the observer are 
among the most variable factors in our sample of BO to B 1 
stars; more so than e.g. T e ff, log(L„), and wind density. 
Inclination angle and i? rot could play an important role 
in determining the precise observational manifestation of 
azimuthally-extended large-scale wind structure that has 
its roots at the photosphere. 

An interesting signature of spiral-shaped structure in 
a hot star wind is the appearance of bowed features. Since 
the winding of a spiral is determined by radial expan- 
sion and stellar rotation, the visibility of bowed features 
in UV time-series data is largely determined by the ratio 
of t'rot/tfoo. Though we only know the projected rotation 
velocity, and our statistical sample is this paper small, the 
ratio of v € sin(i)/t’ co is largest for HD 64740 0.14) and 

HD 157246 (7 Ara; 0.41 for 1995 data set). HD 64760 of 
course exhibits the clearest examples of phase bowing in 
our UV sample (Fig. 10), and indeed it is the best known 
case for this behaviour out of all the UV time-series data 
available for OB stars. Fast rotation is also undoubtedly 
the key to 7 Ara's extreme behaviour, i.e. the transient 
two-component (polar and equatorial) wind. Interestingly, 
HD 150168 is the only other star in our study that shows 
bowed structures which are similar to those evident in 
HD 64760 (Fig. 12). The projected rotation velocity of 
this star is not exceptional and v e sin(t.)/v^. is — 0.09. 
This is possibly then an example then where either the 
inclination angle and/or the density contrast of gas in the 
interacting spiral streams favours the detection of phase 
bowing. 

We have discussed in this paper the best UV time- 
series data sets of early B supergiants available in any of 
the public archives. With the demise of IUE, these col- 
lections will not be improved for the forseeable future. 
One of our key objectives was to collect these homoge- 
neous data sets into a single comparative study. Over 
the next few years these results will be complemented 
by extended time-series studies of luminous early-type 
stars in the Magellanic Clouds using the Far Ultraviolet 
Spectroscopic Explorer satellite (FUSE: AA900 to 1200). 
We anticipate in particular more sensitive searches for ion- 


ization variability in the far-UV lines, and an examina- 
tion of the dependence of wind variability on metallicitv. 
Additionally. —2 metre class optical telescopes are cur- 
rently being used to monitor time-dependent behaviour 
over rotational time-scales in Balmer and He lines of OB 
supergiants. A key goal is these studies is to exploit the 
recombination-dominated lines to determine the coher- 
ence of large-scaie wine structure in the deeper, near- 
star wind regions, clos to potential photospheric inbo- 
rn ogenei ties. 
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